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Abstract 

We compare the observed splitting in the PL spectrum of a strongly coupled 
light-matter system, with the splitting of its dressed modes. In the presence 
of non-negligible decoherence, the two may differ considerably. Whereas the 
dressed mode splitting has a simple expression, the observed splitting has no 
general analytical expression in terms of radicals of the system parameters. 
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There is now a large number of reports of strong light-matter coupling 
in zero dimensional semiconductor systems (one quantum dot in a micro- 
cavity). The successful theoretical description of some of these experiments 
relies on properly taking into account dissipation and decoherence [1]. In 
these approaches, light-matter coupling is described with the linear model 
Hamiltonian H = Uaa'^a + — A)b'^b + g{a^b + ab'^) where a and b are the 
cavity photon and material excitation (bosonic) field operators, respectively, 
with bare mode energies uUa and cj^ — A, coupled with strength g. The pho- 
toluminescence (PL) spectrum of the cavity emission, analyzed in the steady 
state (SS), is greatly aflFected by decoherence. On the one hand, the system 
loses photons and matter excitations at rates 7^ and 75, respectively. On the 
other hand, it is driven by a cw off-resonant pumping or, equivalent ly, by a 
continuous electronic injection in the wetting layer. This is effectively mod- 
elled by two excitation rates. Pa and P5. The dynamics of the system is given 
by a standard master equation. In Ref. 0], we investigated in detail the SS 
dynamics and PL spectrum of this system and obtained general analytical 
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expressions for the lineshapes: 

1 r i-W i + W ^ 
S{uj) = — + = ^{A{u;)} , (1) 

where A{(jj) is a complex function of the real frequency cj, with: 

Q± = uUa — — iV^ ± i? G C (complex frequencies) , 

W =[T_+ i{A/2 + gD)]/R G C (dimensionless) , 

r± = ( r^ ± Tb)/4: G R, Ta^t = ja,b - Pa,b ^ ^ (dccoherence rates) , 

R = ^g'^ -{V_ + i/\/2f G C (half-Rabi frequency) , 

^ _ (at6>^ _ f(7an-7.Pa)(^r+-A/2) Mimensiouless) 
^ - {ata>ss - g^v+{Pa+Pb)+PaVi,{v\+{A/2Y) ^ ^ (.aimensiomessj . 

With these definitions, the strong coupling (SC) regime is achieved when 
the condition |r_| < 5^ is satisfied. In such a regime, the so-called dressed 
modes (DM) appear at the eigenfrequencies overtaking bare modes 

(of the weak coupling (WC) regime) at the frequencies uja and — A. 

The observed spectrum, on the other hand, results from contributions 
from the leading modes (bare in WC or dressed in SC), so its splitting nat- 
urally does not correspond to that of the DM. Emission of each mode is not 
Lorentzian but also has a dispersive part, that can result in an increase or 
decrease of the apparent splitting, or even in its appearance in weak coupling 
or disappearance in strong coupling. It is therefore crucial, experimentally, to 
have an expression for the observed splitting that can be directly compared 
with the observed data. 

It is possible to fit the lineshapes with the expression ([T]) [1]. In this 
text, we focus on the alternative practice of fitting the maxima of these 
spectra, that yield the characteristic crossing or anticrossing behaviors. The 
advantages are in the simplicity (two lines rather than a whole series of 
curves) and the robustness against noise (a parasitic second dot shifts only 
slightly the maxima of the strongly-coupled system whereas its infiuence on 
the whole lineshape can be much more deleterious for the fitting, calling for 
outliers or similar techniques). We need, however, to consider the observed 
splitting, not the dressed mode splitting. We shall see however that this 
approach is hindered by the Abel-Ruffini theorem that prohibits solutions in 
terms of radicals of the parameters. 

For this purpose, we consider the equation dS/duj = 0, that gives the 
local extrema of the spectrum. There exists either one or three real solu- 
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Figure 1: Dressed modes (blue) and observed PL splitting (red) as a function of detuning. 
We fix 76 = 0.1^ and Pa = 0.5g. (a) 7a = 3.8^ and Pb = g, (b) 7a = 3.8^ and Pb = 0.5^, 
(c) 7a = 3.8^ and P5 = 0.1^, and (d) 7a = 4.8^ and P5 = 0.1^. 



tions, corresponding to the singlet or doublet, respectively. Using dS/duj = 
^{dA/duj}^ that holds thanks to G M, we obtain the factorization: 

dSjcu) _ 1 ^{{u; - uj_){uj - co^Xcj - ni)\u; - 1^;)^ _ . . 

duj n ' ^ ^ 

with6j±= b+ -n_+ iW{n^ + n^) ± V(l + (^+ - /2 ^ C. De- 

spite the simple form of Eq. ([2j), none of the roots {uJ±^ f^^} are solutions of it, 
since they are not real. All the solutions must be found taking the imaginary 
part in Eq. ([2j). This yields a quintic equation. It is well known that a quintic 
cannot be solved in general in terms of radicals of its coefficients. This is 
the case of Eq. ([2j) for a general set of parameters. The solutions can still be 
presented in terms of special functions, such as the Jacobi theta functions, 
but such elaborate methods do not serve our objective of providing formulas 
that can be directly used to analyze the experimental data. Therefore, in 
the most general case, we prefer to solve numerically Eq. (ED, exploiting the 
simplicity of the equation rather than the complexity of its solutions. 

In Fig. [U we plot, as a function of detuning: in solid black, the bare 
modes; in dashed black, the DM in the absence of pump and decay; in 
dashed blue, the system DM (3?{f^±}); and in solid red, the maxima in the 
spectra given by numerical solutions to Eq. ([2j). By comparing the red and 
the blue lines, it is clear that both the repulsion of PL lines (cases c and d) 
and their crossing (cases a and b) can appear both in SC (b and c) or WC 
(a and d). 

At resonance, where the spectrum is symmetric with respect to the cen- 
tral frequency cj^, Eq. ([2]) can be solved in terms of radicals. uJa is always 
a solution, corresponding to the minimum in the doublet case and to the 
maximum in the singlet case. The other two possible real solutions give the 
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expression for the observed splitting, given by, in both SC or WC: 




This should be contrasted with the expression for the DM splitting: 

Ac^DM = 2^{R{A = 0)} = 2^3?|Vl-(r-/^)2} . (4) 

Our main statement in this text is that the counterpart of Eq. ([3]) at nonzero 
detuning does not exist in this form and should be computed numerically. 
This is in stark contrast with the DM case, Eq. (jlj), which always has an 
equally simple expression, that in general does not, however, describe accu- 
rately experimental data. 

Expression ([3]) can be used to predict the pumping rates at which the 
splitting will be more visible for a given configuration. The counterpart ex- 
pression for the splitting observed in the direct exciton emission can be ob- 
tained by simply exchanging indexes a ^ b. It is also interesting to note that, 
for the integrated spontaneous PL emission (instead of the SS emission), the 
formula for the observed splitting is found by only removing Pa^b from the F's 
and substituting Pb/Pa by the ratio of initial state populations {b^ b)^ / {a'^ a)^ 
(under the assumption that initially (a'^b)^ = 0). 

In conclusion, we have contrasted the crossing and anticrossing of the PL 
spectrum lines as detuning is varied, with the actual behavior of the dressed 
modes, which truly quantify the strength of the light-matter coupling. We 
showed that a careful analysis is required to properly assess an experiment 
that is not in very strong coupling. No analytical formula in radicals of the 
parameters of the system exists in the general case, supporting the global 
fitting of lineshapes as the most convenient description of experimental data. 
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